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ABSTRACT 

Context. Stellar encounters potentially affect the evolution of the protoplanetary discs in the Orion Nebula Cluster (ONC). However, 
the role of encounters in other cluster environments is less known. 

Aims. We investigate the effect of the encounter-induced disc-mass loss in different cluster environments. 

Methods. Starting from an ONC-like cluster we vary the cluster size and density to determine the correlation of collision time 
scale and disc-mass loss. We use the nbody6++ code to model the dynamics of these clusters and analyze the disc-mass loss due to 
encounters. 

Results. We find that the encounter rate depends strongly on the cluster density but remains rather unaffected by the size of the stellar 
population. This dependency translates directly into the effect on the encounter-induced disc-mass loss. The essential outcome of 
the simulations are: i) Even in clusters four times sparser than the ONC the effect of encounters is still apparent, ii) The density of 
the ONC itself marks a threshold: in less dense and less massive clusters it is the massive stars that dominate the encounter-induced 
disc-mass loss whereas in denser and more massive clusters the low-mass stars play the major role for the disc mass removal. 
Conclusions. It seems that in the central regions of young dense star clusters - the common sites of star formation - stellar encounters 
do affect the evolution of the protoplanetary discs. With higher cluster density low-mass stars become more heavily involved in 
this process. This finding allows for the extrapolation towards extreme stellar systems: in case of the Arches cluster one would expect 
stellar encounters to destroy the discs of most of the low- and high-mass stars in several hundred thousand years, whereas intermediate 
mass stars are able to retain to some extant their discs even under these harsh environmental conditions. 

Key words, stellar dynamics - methods: N-body simulations - stars: pre-main sequence, circumstellar matter 



1. Introduction 

To current knowledge, planetary systems form from the accre- 
tion discs around young stars. These young stars are in most 
cases not isolated but are part of a cluster (e.g. |Lada & Lada 
|2003| . Densities in these cluster environments vary consider- 
ably, spanning a range of 10 pc~ 3 (e.g. rj Chameleontis) to 
10 6 pc 3 (e.g. Arches Cluster). Though it is known that discs dis- 

~||200T] |Hillenbrand|[2002l |Sicilia- 



perse over time (Haisch et al. 



Aguil ar et al.|2 006 1 Curr ie et al.|2008[ ) and that in dense clusters 
(n > 10 J pc" J ) the disc frequency seems to be lower in the core 
(e.g. Balog et al.|2007| l, it is an open question as to how far in- 
teractions with the surrounding stars influence planet formation 
in clusters of different densities. 

Two major mechanisms are potentially able to strongly af- 
fect the evolution of protoplanetary discs and planets in a cluster 
environment: photoevaporation and gravitational interactions. 
Photoevaporation causes the heating and evaporation of disc 
matter by the intense UV radiation from massive stars. First 
model s of p hotoevaporation were develope d by |Johnstone et al.| 
(1998| and |St6rzer & Hollenbach| ( fT999| l (see also references 
in Hollenbach, Yorke & Johnstone 2000), and have been much 
improved in the past years (Cla rke et al.||2001| |M atsuyama 
et al.|20031|Joh nstone et al.|2004HAlexander et al'|2005| |2 006 ; 



Ercol ano et al.||20081 |Gorti & Hollenbach||20091 [Drake et al.| 
2009). Gravitational interactions are another important effect on 
the population of stars, discs, and planets in a cluster environ- 
ment. An encounter between a circumstellar disc and a nearby 
passing star can lead to significant loss of mass and angular mo- 



mentum from the disc. While such isolated encounters have been 
studied in a large variety (|Heller|19"9"3"j|Clarke & Pringle|1993| 
I Qstrikerl 1 9941 |Heller| 1995) |Hall et al. 1 1 9961 |Hall| 19971 IPfalzner 
20041 |Pfalzner et al.|2005t|Moeckel & Bally|2006||2007bl|Kley 
|et al.|2008 1, only few numerical studies have investigated the ef- 
fect of stellar encounters on circumstellar discs in a dense cluster 
environment directly ( Scally & Clarke 2001 Adams et al. 2006). 

Only recently has it been shown from numerical simulations 
that stellar encounters do have an effect on the discs surrounding 



stars in a young dense cluster ( Olczak et al.| 2006 ; Pfalzner et al 
2006 ; |Pfalzner 2006; M oeckel & Bally|2006| |2007a|b| |Pfalzner 



& Olczak 200 7a|b| see also the review by Zinnecker & Yorke 
2007). The massive stars in the centre of such a stellar clus 



ter act as gravitational foci for the lower mass stars (Pfalzner 
let al.||2006b. Discs are most affected when the masses of the 



stars involved in an encounter are unequal ( Olczak et al.| 2006 



|Moeckel & Ba lly 2007b), so it is the massive stars that domi- 
nate the encounter-induced disc-mass loss in young dense clus- 
ters ([Olczak et al. 2006). Observational evidence for this effect 



Olczak et al. (2008). 



has been found by 

The numerical results obtained in our previous investigations 
are based on a dynamical model of the Orion Nebula Cluster 
(ONC) - one of the observationally most intensively studied 
young star cluster. It was demonstrated that in the ONC stel- 
lar encounters can have a significant impact on the evolution 
of the young stars and their surrounding d iscs (j Pfalzner 2006; 
[ Olczak et al.||2006l |Pfalzner et al.|[2006l |Pfalzner & Qlczak| 
|2007a||Olczak et al.|2008HPfalzner et al.|2008] >. However, inves- 
tigating one model star cluster is not sufficient to draw general 
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conclusions - in fact, one could not answer questions as: How 
would things change in a denser cluster? Would a higher density 
inevitably imply that stellar encounters play a more important 
role in the star and planet formation process? And what would 
be the situation in more massive clusters? Would the larger num- 
ber of stars - in particular massive stars - play a role? 

Conclusive answers to these questions demand further nu- 
merical investigations covering a larger parameter space in clus- 
ter parameters. Here we realise this by modelling scaled ver- 
sions of the standard ONC model - clusters with varying stel- 
lar numbers and sizes. The focus of this investigation is on the 
encounter-induced disc-mass loss. Throughout this work we as- 
sume that initially all stars are surrounded by protoplanetary 
discs. This is justified by observations that reveal disc fractions 
of nearly 100% in very young star cluster s (e.g. |Haisch et al. 
[20001 |Lada et^|2000l |Haisch et al.|[200T] |Hillenbrand||2005) . 
In Section |2J we briefly outline the observationally determined 
basic properties of the ONC, that serves as a reference model 
for the other cluster models. The computational method is de- 
scribed in Section [5] and the construction of the numerical mod- 
els is outlined in Section [4] Afterwards we present results from 
a numerical approach to this problem in Section [5] and compare 
with analytical estimates in Section [6] The conclusion and dis- 
cussion mark the last section of this paper. 



2. Structure and dynamics of the ONC 

The ONC is a rich stellar cluster with about 4000 members 
with masses m > 0.08 M Q and a radial extension of ~2.5 pc 



( |Hillenbrand & HartmamifT99 8 ; Hillenbran d~& Carpenter|2 000 ) 
Most of the objects are T Tauri stars. The mean stellar 
mass is about 0.5 Mo and the half-mass radius Rh m ~ 1 pc 
( |Hillenbrand & Hartmann||1998|l. Recent studies of the stellar 
mass distribution (Hillenbrand & Carpenter 2000, Luhman et al. 
2000} |Muench et al.|2002||Slesnick et al.|2004|> reveal no si gnif- 
icant deviation from the generalized IMF of |Kroupa| ( |2002[ ), 



£(m) 



mr 13 , 0.08 < m/Mo < 0.50, 
0.50 < m/M < 1.00, 



-2.3 
-2.3 



(1) 



1.00 < m/M < oo 



The mean age of the whole cluster has been estimated to be 
7onc ~ 1 Myr, although a significant age spread of the individual 
stars is evident ( [Hillenbrand[l997l|Palla & Stahler|2000) . 

The density and velocity distribution of the ONC resembles 
an isothermal sphere. The central n umber density p mre in the 



inner 0.053 pc reaches 4.7 ■ 10 4 pc" 3 ( jMcCaughrean et al.|2(K)2 1 



and makes the ONC the densest nearby (< 1 kpc) young stel- 
lar cluster. The dense inner part of the ONC, also known as the 
Trapezium Cluster (TC), is characterized by 7?xcS0.3pc and 
N TC ~ 500, or pxc ~ 4 ■ 10 3 pc" 3 . 

In the most detailed study on circumstellar discs in the 



Trapezium Cluster, Lada et al. (2000) found a fraction of 80- 
85 % discs among the stellar population from the L-band excess. 
This is in agreement with an earlier investigation of the complete 
ONC in which [Hillenbrand et al. ( 1998 ) report a disc fraction of 
50-90 % (though relying only on Iq - K colors) and justifies the 
assumption of a 100% primordial disc fraction in the here pre- 
sented simulations. 

In the following we describe the construction of the numeri- 
cal cluster models that have been used in our simulations. 



3. Computational method 

The basic dynamical model of the ONC used here is described 
in |Qlczak et al.| (|2006|), w ith several extensions discussed in 
Pfalzner & Olczak| ( |2007b[ ). We summarize the main aspects of 
our model: The initial stellar population consists of 4000 mem- 
bers with masses between 0.08 M and 50 M sampled from the 
standard Kroupa IMF (see Eq. (|T]i). The system is initially in 
virial equilibrium, 



2GM 



0.5. 



(2) 



where Rh m is the half-mass radius of the cluster, <x the velocity 
dispersion, M the total mass, and G the gravitational constant. It 
is characterized by a radial density profile, p cc r~ 2 , with a central 
density ~4T0 4 pc~ 3 . We adopt a Maxwellian velocity distribu- 
tion as would be expected from theory of star cluster formation 
(e.g. Clar ke et al.p OOO) and roughly in agreement with obser- 
vations of star-forming regions (e.g. Falgarone & Phillips|1990 



Mies ch & Scalo|1995| l. We use here a single star model only and 
do not include the effect of gas expulsion or stellar evolution. All 
simulations have been performed with the direct Af-body code 
nbody6++ ( |Spurzem| 1 999) |Aarseth|2003) . 

For the generation of star cluster models in the present in- 
vestigation the initial radial density profile has been modified 
slightly. To first order the isothermal sphere, p(r) cc r~ 2 , repre- 
sents the projected density distribution of the ONC, yet a flatten- 
ing in the core, E(r) oc r~ 5 , < r < R cme , R cole ~ 0.2 pc, is ob- 
served ( Scally et al.|2005 1. Validating the initial setup by means 
of the best reproduction of the current projected density distribu- 
tion of the ONC after a simulation time of 1 Myr, the evaluation 
of numerous initial configurations led to the following best esti- 
mate of the initial three-dimensional density distribution: 



Po(r) 



Po(r/R cme r 23 



P0 (r/Rcom) 





-2.0 



r G (0,/? CO re] 
f 6 (Rcote,R] 

r e (R, oo] 



(3) 



where po = 3. 1 • 10 3 stars pc~ 3 , R cme = 0.2 pc, and R = 2.5 pc. 

Moreover, the generation of the high-mass end of the mass 
function has been modified. In the case of the ONC the upper 
mass was chosen to be 50 M because this value corresponds 
to the mass of the most massive stellar system in the ONC. 
However, stars with larger masses are expected to form in more 
massive clusters (Oey & Clarke 2005 , Weidner & Kroupa 2006 ). 
Thus in the framework of this numerical investigation the up- 
per mass limit has been set to the current accepted fundamental 
upper mass Umit, m ma x = 150 M (fFT ger 2005 ; |Oey & C larke 
~ Koen|p006l | Weidner & Kroupa||2006| |Maiz Apellaniz 



2005 



etal. 



2007[|Zinnecker & Yorke[2 007 ). 



The choice of a fixed upper mass limit, though in disagree- 
ment with the well-established non-trivial correlation of the 
mass of a star cluster and its most massive m ember (e.g. 
2003| |Weidner & Kroupap)06l paschberger & Clarke 



Larson 
2008), 



known. 



was motivated by the fact that the exact relation is not 
However, a comparison with the "sorted sampling algorithm" of 
Weidner & Kroupa (2006) in Table [T] shows that - at least in a 
statistical sense - the exact prescription for the generation of the 
maximum stellar mass in a cluster is not as important as it might 
seem. The values obtained by random sampling are only slightly 
larger than those from sorted sampling. 

Stellar encounters in dense clusters can lead to significant 
transport of mass and angular momentum in protoplanetary discs 
dOlczak et aTJ|2006] [Pfalzner et al.|[2006l |Pfalzner & Qlczak| 
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family of models 
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[pc] 


[10 3 pc" 3 ] 


[km s-'] 




[pc] 


[10 3 pc" 3 ] 


[kms" 1 ] 


[Mo] 


[M ] 


[M G ] 


1000 


DO 


2.50 


1.3 


1.15 


SO 


0.63 


4.8 


2.42 


36 


32 


25 ± 


15 


2000 


Dl 


2.50 


2.7 


1.64 


SI 


1.25 


5.1 


2.37 


52 


45 


25 ± 


15 


4000 


D2 


2.50 


5.3 


2.26 


S2 


2.50 


5.3 


2.26 


79 


63 


55 + 


25 


8000 


D3 


2.50 


10.5 


3.11 


S3 


5.00 


5.3 


2.13 


94 


80 


75 ± 


25 


16000 


D4 


2.50 


21.1 


4.34 


S4 


10.00 


5.3 


2.11 


125 


112 


95 ± 


35 


32000 


D5 


2.50 


42.1 


6.04 


S5 


20.00 


5.2 


2.03 


137 


126 


100± 


35 



Table 1. Averaged initial parameters of the cluster models, divided among the families of density-scaled and size-scaled models. 
The first column contains the number of stars, N, the next eight columns contain the designation, the cluster size, R, the number 
density in a sphere of radius 0.3 pc, p-rc (equivalent to the Trapezium Cluster in the ONC, see Section[2]), and the three-dimensional 
velocity dispersion, cr 3D , for the density-scaled and size-scaled models, respectively. The last three columns denote the median 
maximum stellar mass in each simulation, m™^, t he median maximum stella r mass for sorted sampling, m™^ x , and estimates from 



observational data, f7j^ a s x , both taken from Fig. 7 of Weidner & Kroupa (2006 1. 



2007a). In the present investigation we assume that all discs 
are of low mass, i.e. a low mass ratio of disc and central star, 
m disc/ m * 0.1, in agreement with observations of the ONC 



( [Bally etal.|199 8 ; Willia ms et aL]2 005 ; Mann & Williams 2009) 
We use Eq. (1) from Pfalzner et al. '( 2006) > to keep track of the 
relative disc-mass loss of each star-disc system due to encoun- 
ters. Approaches of stars are only considered to be encounters if 
the calculated relative disk-mass loss is higher than 3 %, corre- 
sponding to the lcr error in our simulations of star-disc encoun- 
ters. In this case the relative disc-mass loss is independent of the 
disc mass and depends only on the mass ratio of the interact- 
ing stars and the orbital parameters (see [Pfalzner et alT 2006). 
Our estimate of the accumulated disc-mass loss is an upper limit 
because the underlying formula is only valid for co-planar, pro- 
grade encounters, which are the most perturbing. A simplified 
prescription assigns stars into one of two distinct groups: if the 
relative disc-mass loss exceeds 90 % of the initial disc mass, 
stars are marked as "discless"; otherwise they are termed "star- 
disc systems". As in our previous investigations, the determi- 
nation of the disc-mass loss involves two different models of 
the initial distribution of disc sizes: i) scaled disc sizes with 
r& = 150 AU yjM\ [M Q ], which is equivalent to the assumption 
of a fixed force at the disc boundary, and ii) equal disc sizes with 
= 150 AU. Whenever results are presented, we will specify 
which of these two distributions has been used. 



4. Numerical models 

For the present study we have set up a variety of scaled ver- 
sions of the standard ONC model cluster with varying sizes and 
stellar numbers. In total, eleven cluster models have been set up 
(see Table [TJ. They form two parametric groups, the "density- 
scaled" (D0-D5) and the "size-scaled" (S0-S5) group, both con- 
taining six clusters with stellar numbers of 1000, 2000, 4000, 
8000, 16000, and 32000, respectively. Models D2 and S2 are 
identical and correspond to the standard ONC model with the 
adopted higher stellar upper mass limit, m max = 150 M Q . The 
other ten cluster models have been set up as scaled represen- 
tations of ONC-like clusters. As in the case of the numerical 
model of the ONC, for each cluster model a set of simulations 
has been performed with varying random configurations of po- 
sitions, velocities, and masses, according to the given distribu- 
tions, to lower the effect of statistical uncertainties. For the clus- 
ters with 1000, 2000, 4000, 8000, 16000, and 32000 particles, 
a number of 200, 100, 100, 50, 20, and 20 simulations seemed 
appropriate to provide sufficiently robust results. 



It has to be noted that these "artificial" stellar systems are 
not just theoretical models but have as well counterparts in the 
observational catalogues of star clusters: the young star cluster 
NGC 2024 (e.g. |Haisch et al.pOOOl |Bik et al.pOOlj |Liu et al. 
2003) is well represented by the 1000 particle model, whereas 
the 16000 particle model has its counterpart in the massive clus- 
ter NGC 3603 (e.g. |de Pree et al.|1999[[Stolte et al.|2004[[2006 1. 



Density-scaled cluster models 

The six density-scaled cluster models (D0-D5) have been sim- 
ulated with the same initial size as the ONC (R = 2.5 pc). Due 
to the adopted number density distribution, roughly represented 
by p(r) = po r~ 2 , the density of the models scales as the stellar 
number, 



N 



f 

Jo 



p(r) r dr d£l oc p R , 



(4) 



though for an exact treatment one would have to consider the 
steeper density profile of the core, pcoreM = Pcore.o r 



-2.3 



-f 

Jo 



Pcore(r) r drdtlcc p core R 



,0.7 



(5) 



However, since the core population is not dominant in terms of 
number, the clusters are characterised in good approximation by 
densities that are 1/4, 1/2, 1, 2, 4, and 8 times the density of 
the ONC (at any radius), respectively. These models are used 
to study the importance of the density for the effect of star-disc 
encounters in a cluster environment. 



Size-scaled cluster models 

Five more cluster models have been simulated with the same 
initial density as the ONC but with varying extension. The set 
of size-scaled cluster models (S0-S5) is used to study the pure 
effect of the size of the stellar population. Due to the relation Q, 
the initial size of these clusters scales as the stellar number and 
was set up with 1/4, 1/2, 1, 2, 4 and 8 times the initial size of the 
ONC, respectively. 

The initial parameters of the cluster models, for each model 
averaged over all configurations, are presented in Table [T] Here 
the number density in the Trapezium Cluster, p T c, is taken as 
a reference value for all simulations. As expected, the density 
scales with the number of stars for the density-scaled models, 
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while it is rather constant for the size-scaled models. The veloc- 
ity dispersion, that satisfies the relation 



2GM 



R 



(6) 



shows the expected scaling of V^V for the density-scaled models, 
and is again roughly constant for the size-scaled models, as ex- 
pected from N cc R (Eq. (|4 1) and the above relation. The reason 
for the slight increase of the velocity dispersion with decreas- 
ing stellar number for the size-scaled models is the steeper den- 
sity profile in the cluster core, which becomes more dominant in 
terms of stellar number with decreasing cluster size. Combining 
N core °c R 01 from Eq. ^ and the above relation gives roughly 
Ccore K N~ 03 and thus explains the correlation. 
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5. Results of numerical simulations 

In this section the results of the numerical simulations of the two 
families of cluster models will be presented. A short discussion 
of the characteristic scaling relations and differences of the clus- 
ter dynamics will be followed by a more detailed investigation 
of the encounter-induced disc-mass loss among the two families 
of density-scaled and size-scaled cluster models. 



Fig. 1. Projected density profiles of the density-scaled models 
compared to observational data. The initial profile (grey lines) 
and the profile at a simulation time of 1 Myr (black lines) are 
shown. From bottom to top in each colour regime the cluster 
models D0-D5 are marked by a short-dashed, long-dashed, solid, 
dotted, dot-long-dashed, and dot-short-dashed line, respectively. 
The observational data are from a compilation of McCaughrean 
[eTaT1 ( |2002l i and |mienbrandl ( fl997) . 



5.1. Cluster dynamics 

In Fig. [T]the evolution of the projected density distribution of 
the density-scaled models is shown. The shape of the distribu- 
tions is in all cases very similar. The evolved distributions have 
nearly identical shapes and are separated by vertical intervals 
of 0.3 in log-space, which corresponds to the difference of the 
initial densities by a factor 2. Only in the innermost cluster re- 
gions slight deviations between the evolved distributions are ap- 
parent. These are attributed to the poorer random sampling of the 
initial particle distribution due to the very steep density profile, 
p(r) oc r~ 2 - 3 , as is evident from the larger scatter among the ini- 
tial profiles. However, after 1 Myr these deviations are smoothed 
out to a large degree. 

Due to the nearly exact qualitative and scaled quantitative 
evolution of the density-scaled cluster models, it is justified to 
ascribe differences of the effects of star-disc encounters on the 
stellar population mainly to one parameter, namely the initial 
density of the cluster models. Nevertheless, the influence of the 
particle number has to be considered, too. 

The size-scaled cluster models show a different dynamical 
evolution compared to the density-scaled models. The tempo- 
ral evolution of the densities pxc in Fig. [^demonstrates that the 
clusters evolve on slightly different time scales, where the den- 
sity declines faster for the less populated clusters SO and SI. 
However, the densities of the models S1-S5 differ not much, and 
are consistent with a coeval decline. The evolution of the clus- 
ter densities does not - to first order - depend on the number 
of particles, probably with the exception of the model SO. The 
size-scaled models are thus well suited to investigate the effect 
of the number of cluster stars on star-disc encounters and the 
corresponding induced disc-mass loss. 

5.2. Encounter dynamics 

We investigate the effect of the cluster density on the encounter- 
induced disc-mass loss via the evolution of the cluster disc frac- 
tion (CDF). The distributions in Fig. [3] show the average frac- 




1 2 3 4 5 



simulation time [Myr] 

Fig. 2. Time evolution of stellar densities pjc (in a volume of 
radius Rjc = 0.3 pc, see Section [2| of the size-scaled clus- 
ter models. The cluster models S0-S5 are marked by a short- 
dashed, long-dashed, solid, dotted, dot-long-dashed, and dot- 
short-dashed line, respectively. The horizontal error bar marks 
the corresponding observational estimates for comparison. 



tion of stars that are surrounded by disc material and corre- 
spond from top to bottom to clusters with increasing density. 
The curves have been smoothed with Bezier curves to provide 
a clearer view. It is evident that the fraction of destroyed discs 
in the Trapezium Cluster increases significantly with increasing 
cluster density. In particular, the effect becomes much stronger 
for the models D3-D5 with 2-8 times the density of the ONC. 
In the case of the densest model D5, even up to 60% of the 
stars in the Trapezium Cluster could have lost their discs after 
1 Myr of dynamical evolution. But it is also interesting to note 
that even in a cluster 4 times less dense than the ONC (model 
DO), still 10-15 % of the stars could lose their surrounding discs 
due to gravitational interactions with cluster members. However, 
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0.3 



2 3 
simulation time [Myr] 



Fig. 3. Time evolution of the cluster disc fraction of the density- 
scaled models for a region of the size of the Trapezium Cluster 
(Rtc = 0.3 pc). The curves have been smoothed by Bezier curves 
to avoid intersecting lines. Here initially equal disc sizes have 
been assumed. From top to bottom the cluster models D0-D5 are 
marked by a short-dashed, long-dashed, solid, dotted, dot-long- 
dashed, and dot-short-dashed line, respectively. The typical error 
bar is indicated in the upper right. 



one has to treat these numbers with care due to the - partially 
significant uncertainties that go into the calculations (see |Qlczak 
et al.|2006| ). Nonetheless, what is more important here - and re- 
lies only on the relative quantities - is the fact that the distribu- 
tions in Fig.[3]are not equidistant but do show larger differences 
with increasing density. 

The impression that the disc fraction increases at later stages 
of the cluster evolution is solely due to the acceleration of per- 
turbed systems leaving the sampling volume. The escape rate 
follows the trend of the disc destruction rate with a time delay 
that is determined by the difference of the crossing time of disc- 
less stars and star-disc systems in the sampling volume. Hence, 
the decrease of the disc fraction is followed by an increase and 
eventually remains constant some time after disc destruction has 
stopped. 

The preferred mass range of stars that become dispropor- 
tionately high involved in perturbing encounters changes signif- 
icantly between the models D2 and D4. From previous inves- 
tigations (Olcz ak et al.| 2006 Pfalzner et al.||2006) l we already 
know that in the model D2, representing an ONC-like cluster, 
the high-mass stars' discs are mostly affected by encounters. 
These encounters occur preferentially with low-mass stars due 
to gravitational focusing. However, in the model cluster D4, a 
four times denser system, it is the low-mass stars that dominate 
the disc-mass loss of the cluster population. 

We demonstrate this transition by binning the fraction of en- 
counters, normalised to the total number of encounters of each 
model (so that the integrated surface is unity in each case). 
Recall that encounters have been defined as flybys causing more 
than 3 % disk-mass loss (see Section [3]). Fig. [4] shows for the 
group of low-mass stars that the dominance of perturbations 
from much more massive stars in the models D0-D2, represented 
by the peak at relative perturber masses > 100, changes dramati- 
cally towards encounters of nearly equal-mass stars in the model 
D4 (with a peak at relative perturber masses m 3). Similarly, the 
typical geometry of the encounter orbits underlies a strong tran- 
sition. As is evident from Fig. [5] most encounters of low-mass 
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Fig. 4. Comparison of the fraction of encounters as a function of 
the relative perturber mass, i.e. the mass ratio of perturber and 
perturbed star, of the group of low-mass stars (see Appendix |A| 
for the width of the mass intervals) for the models DO (dotted), 
D2 (solid), and D4 (dot-long-dashed), respectively. Here initially 
equal disc sizes have been assumed. 
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Fig. 5. Comparison of the fraction of encounters as a function 
of eccentricity of the group of low-mass stars (see Appendix |A| 
for the width of the mass intervals) for the models DO (dotted), 
D2 (solid), and D4 (dot-long-dashed), respectively. Here initially 
equal disc sizes have been assumed. 



stars in the models D0-D2 are parabolic (with an eccentricity 
e * 1). In contrast, in the model D4 these stars usually expe- 
rience a strongly hyperbolic flyby of a perturber, mostly with 
eccentricity e > 10. Hence, with increasing cluster density, the 
dominant mode of star-disc interactions changes towards hyper- 
bolic encounters of low-mass stars between each other. 

This transition is a consequence of gravitational focusing by 
high-mass stars becoming significantly less important among the 
models D0-D5. Gravitational focusing forces a deflection of the 
interacting stars onto a less eccentric orbit, effectively increasing 
the cross-section above geometrical: 



where 







(Vrel) 2 



ft V min 

16 o- 2 



(7) 



(8) 
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is the gravitational focusing term or Safronov number, b the im- 
pact parameter, r mm the minimum distance, v m i„ the escape speed 
at the minimum distance, (v re i) the mean relative speed, and <x 
the cluster velocity dispersion (one finds (v re i) = Act I -0r for a 
Maxwellian distribution with dispersion <x). Adopting the typi- 
cal cluster and encounter parameters from Table [T] and [2] respec- 
tively, and substituting v m ; n = v enc and r m ; n = r enc , we find that 
gravitational focusing by low-mass stars is negligible (0 <K 1) 
in all our cluster models. In contrast, high-mass stars in the 
model DO play a substantial role as gravitational foci (0 * 120), 
whereas their cross section is reduced by more than one order of 
magnitude in the model D5 (0 ~ 4). 

In summary, the encounter-induced disc-mass loss in clus- 
ter environments of different densities shows two important fea- 
tures: i) low-mass stars lose a larger fraction of their disc mate- 
rial with increasing cluster density, and ii) the discs of the most 
massive stars are (nearly) completely destroyed, independent of 
the density of the cluster environment. The important finding for 
i) is that the correlation is not linear, but shows a much larger 
increase for the cluster models with higher densities than model 
D2, implying that there exists a critical density p cr it that marks 
the onset of a much more destructive effect of star-disc encoun- 
ters. This critical density seems to be close to the density of the 

ONC, p c ,it ~PONC- 

We find that the evolution of the CDF in the Trapezium 
Cluster region is very similar among the size-scaled models and 
corresponds to the distribution of the model D2 in Fig. [3] The 
size-scaled models are obviously equivalent in their environmen- 
tal effect on protoplanetary discs despite the slightly different 
dynamical evolution. The density of the models SO and SI de- 
creases faster than for the more massive clusters, even up to a 
factor of 2 in case of the model SO (see Section |5T) . Thus one 
would expect a lower encounter rate in these smaller systems 
and, consequently, on average a lower disc-mass loss. However, 
this is obviously not the case, and is explained by the fact that, 
similarly to the finding for the density-scaled models, high-mass 
stars become less important as gravitational foci for the low- 
mass stars in clusters with larger stellar populations. Hence in 
terms of encounter statistics, the lower density of the models SO 
and SI is compensated by the more frequent interactions of the 
high-mass stars. 

5.3. Validation of the numerical method 

So far the disc-mass loss has been calculated from Eq. (1) of 
Pfalzner et al. (2006) by treating all encounters as parabolic. To 
account for the lower disc-mass loss in hyperbolic encounters, 
we have determined a function that quantifies the reduction of 
the disc-mass loss in dependence of the orbital eccentricity: 



Am(e) =exp[0.12(e - 1)] x 
{0.83 -0.015(e- 1) 



0.17exp[0.1(e- 1)]}. 



(9) 



Eq. (j9j is a fit function to the median distribution of the rela- 
tive disc-mass loss as a function of eccentricity, normalised to 
the parabolic case, for all star-disc simulations that have been 
performed. 

As shown in Fig. [6] the number of hyperbolic encounters 
is significantly reduced if the eccentricity is considered explic- 
itly in the calculation of the disc-mass loss using the fit func- 
tion This is a consequence of our definition of an encounter 
(see Section [3J: the fraction of perturbations that cause a rela- 
tive disc-mass loss above 3 % is lower for higher eccentricities, 
hence the number of eccentric encounters decreases. However, 
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Fig. 6. Number of encounters as a function of eccentricity (log- 
arithmic bins), plotted for three different groups by means of 
disc-mass loss per encounter. The white surface represents all 
encounters (i.e. minimum 3 % disc-mass loss, cf. |Qlczak et al.| 
2006 ), the light grey those that removed at least 50 % of the disc- 
mass, while the dark grey stands for the most destructive encoun- 
ters that caused a disc-mass loss of at least 90 %. The two plots 
show the distributions for the model D4. Top: Disc-mass loss cal- 
culated assuming parabolic encounters. Bottom: Disc-mass loss 
calculation corrected for effects of eccentricity by using Eq. ([9]). 



because the fit function represents the median distribution of all 
simulated star-disc encounters and the effect of strongly perturb- 
ing encounters is only weakly dependent on the eccentricity, the 
number of strongly perturbing encounters in Fig. (light and 
dark grey surfaces) is underestimated. 

Hence it is the disc-mass loss induced by weak hyperbolic in- 
teractions that has been overestimated in the calculations. These 
events are most numerous in the models D3-D5 with densities 
P > Pciit and result preferentially from interactions of low-mass 
stars with roughly equal-mass perturbers. Consequently, when 
considering explicitly the reduced disc-mass loss in hyperbolic 
encounters (see Fig. [7]), the outstanding role of these dense clus- 
ters as environments of huge disc destruction becomes less pro- 
nounced (cf. Fig. |3j, though the encounter-induced disc-mass 
loss is still considerably larger compared to the sparser clusters 
D0-D2. 
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Fig. 7. Time evolution of the cluster disc fraction of the density- 
scaled models, not restricted to parabolic encounters, for a re- 
gion of the size of the Trapezium Cluster (Rjc = 0.3 pc). The 
curves have been smoothed by Bezier curves to avoid intersect- 
ing lines. Here initially equal disc sizes have been assumed. 
From top to bottom the cluster models D0-D5 are marked by 
a short-dashed, long-dashed, solid, dotted, dot-long-dashed, and 
dot-short-dashed line, respectively. The typical error bar is indi- 
cated in the upper right. 
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Fig. 8. Average relative disc-mass loss at 1 Myr for the 
Trapezium Cluster as a function of the stellar mass for initially 
equal disc sizes. The standard ONC model with a stellar upper 
mass limit of 50 M Q (dark grey bars) is compared to a simu- 
lation of the same model with an upper mass limit of 150 M 
(light grey bars). 



The effect of the higher upper mass limit on our simula- 
tion results compared to the basic ONC model is investigated 
in Fig. [8] It shows the average relative disc-mass loss as a func- 
tion of the stellar mass for the standard ONC model with a stel- 
lar upper mass limit of 50 M (dark grey boxes), and the same 
model with an upper mass limit of 150 M (light grey boxes). 
For masses below 10 M the two distributions are quantitatively 
in good agreement. 

In the range of ~ 10-50 Mo the average disc-mass loss of the 
50M o -limit model is significantly higher. This is to be expected 
because stars in this mass range can act as additional strong grav- 
itational foci in the presence of a 50 M star, while their effect 



is largely reduced if a 150 M star is gravitationally dominat- 
ing. The most massive stars in the range 80 - 15OM show, as 
expected, the largest average relative disc-mass loss. That it is 
somewhat lower than that of the most massive star of the 50 M - 
limit model is in agreement with the stronger gravitational at- 
traction of their disc, leading on average to a reduced disc-mass 
loss per encounter. However, since the highest mass bin is only 
populated by 9 stars, any further conclusions about the average 
relative disc-mass loss of the most massive stars would be highly 
speculative. 

6. Comparison of numerical results and analytical 
estimates 

In this section the numerical results will be compared to analyt- 
ical estimates. The treatment of encounters involves one impor- 



tant time scale, the collision time f co u (Binney & Tremaine 1987 
Eq. 8-123), 



coll 



16 y/nparl 1 + 



Cm* 



i 



(10) 



2<x 2 r lk 

Here the inverse of the collision time will be introduced as the 

f coir 



encounter rate, / enc = t Using the escape velocity v* from the 



stellar surface, 



v* = 



2Gm+ 



(11) 



the encounter rate can be written as 

/ enc = 16Vip<^l + ^j, (12) 

where G denotes the gravitational constant, m* and r* the stellar 
mass and radius, and p and cr the density and velocity disper- 
sion of the star cluster. In the following, the stellar radius r+ will 
be replaced by the "typical interaction radius" r enc , that means 
the radius at which the star is subject to a significant (but still 
frequent) perturbation that potentially can remove some frac- 
tion of the disc-mass. Eq. ( 12 1 will be evaluated for three dif- 



ferent stellar masses, representing stellar mass groups of low-, 
intermediate- and high-mass stars. Appropriate typical interac- 
tion radii have been taken from Table 3 of Olcz ak et al.| (|2006). 
The set of masses m+, radii r enc , and "encounter escape speeds" 
Venc resulting from Eq. ( 1 1 1 is shown in Table|2] The last column 
contains the "gravitational focusing parameter" y enc , an approx- 
imation parameter defined as the power of ten best representing 



7e 



= ioLl°gv e 2 nc + l/2J 



(13) 



where \_x\ denotes the floor function 1 of x. For clarity, mass- 
dependent quantities, x enc , will be explicitly denoted as functions 
of mass, x enc (m*), in the following. 

Because the model D2/S2 represents the standard ONC 
model, which has been intensively studied, the calculations will 
be normalised to this model. All quantities related to this model 
will be thus denoted by a "0" as subscript. Adopting the ini- 
tial velocity dispersion of the model D2/S2, <x as 2.3kms _1 w 
0.5 AU/yr, using 



4o- z = Ao-i 



AU 2 
yr 2 



The floor function LxJ gives the largest integer less than or equal to 
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mass groups 


m 


* 


' enc 


Vcnc 


Tenc 




[M 


g] 


[AU] 


[AU/yr] 




low mass 


0.1 - 


1 


1CM 


-0.3 


10-' 


intermediate mass 


1 - 


10 


10 2-2.5 


~ 2 


1 


high mass 


10- 


100 


10 2 


~ 6 


10 



Table 2. Typical parameters adopted for the calculation of the 
encounter rate of cluster stars. The first column denotes the 
three mass groups, the second column contains the adopted mass 
ranges, m*, while in the third typical interaction radii, r enc , are 
listed for each mass group. In the last two columns the resulting 
"encounter escape speeds", v enc , and the "gravitational focusing 
parameter", y enc , are noted (see text for definitions of these quan- 
tities). 



and the numbers given in Table [2j Eq. ( 12 1 can be simplified to 



/enc = 16 V^PO-' - enc( m *) 



1 +yenc('«*)(^) 



(14) 



An even more compact representation is achieved by consider- 
ing the scaling properties of the two families of models: as can 
be derived from Eqs. Q and |6|, the scaling relations for the 
density-scaled models are p oc N and <x oc yfN, while p = const 
and <x = const is found for the size-scaled models. Using these 
relations, transforming r enc (m*) to y e nc('w*) via Eqs. (|TTJ and 
( pj) , and normalising the encounter rate to the model D2/S2, 

/enc™ = /enc//enc,0, One obtains 



f 
■u 



norm 
enc 



m* y e nc('M*,o) 

m*,0 7enc('«*) 

PQ J 1 + Tenc('«*,0) 

1 + 7enc(™*) 



1 + renc(w* >0 ) 



D-models , 



S-models . 



(15) 



The derived relation for the normalised encounter rate pre- 
dicts very different scaling relations for the two families of clus- 
ter models. Density scaled-models are expected to show large 
variations of the number of encounters, with a strong depen- 
dency on the density for low-mass stars, i.e. when y e nc('«*) ^ 1. 
In contrast, the encounter rate for the size-scaled models is ex- 
pected to vary only for the different mass groups, but not among 
different models. For a better overview of the scaling in terms 
of numbers, Table [3] lists approximated relative encounter rates 
/enc™' normalised to the low-mass group of the model D2/S2. 
Table [3] demonstrates that the gravitational focusing parameter 
plays an important role for the massive stars. The encounter rates 
increase dramatically by roughly one order of magnitude from 
the low- and intermediate-mass stars to the high-mass stars for 
the model D2/S2. This finding agrees well with the number of 
encounters of the ONC model found earlier (cf. Pfa lzner et al.| 
2006) . 

In summary, what one would expect from the numerical sim- 
ulations of the density-scaled models is a steep increase of the 
encounter rate oc p 3 ' 2 in case of the low-mass stars and a consid- 
erably shallower dependency oc p l l 2 for the high-mass stars. For 
low densities, corresponding to low particle numbers, one would 
expect that the high-mass stars dominate the encounter rate via 
gravitational focusing, favouring on overall scaling oc p l l 2 . In 
contrast, the size-scaled models should produce very similar re- 
sults in terms of encounter rate, independent of the specific clus- 
ter model. 



These expectations are in good agreement with the results 
from the numerical simulations. We demonstrate this in Fig. [9] 
via the average encounter rate of stars of all masses and of the 
three mass groups, respectively, normalised in each case to the 
model D2/S2. Fig. [9^ shows that the encounter rates of the clus- 
ter models D0-D2 are scaling roughly as N 1 ^ 2 . For higher par- 
ticle numbers the distribution becomes more complex. Here the 
high-mass stars show a trend of decreasing encounter rate with 
particle number. This feature accounts for the decreasing impor- 
tance of the high-mass stars as gravitational foci (for the lower 
mass stars) and is a consequence of the decreasing ratio of the 
mass of the most massive star and the cluster mass. Accordingly, 
the distribution of the encounter rate tends towards the analytical 
limit of N 3 ^ 2 for low-mass stars, representing the more frequent 
interaction of low-mass stars with each other in the models D3- 
D5. 

Fig. [9J3 shows that the size-scaled models are equivalent in 
their environmental effect on protoplanetary discs. In the case of 
low- and intermediate mass stars the presented encounter rates, 
normalised to the model S2, are in good agreement with the ana- 
lytical estimate, which predicts a constant distribution as a func- 
tion of particle number. In contrast, the normalised encounter 
rate of the high-mass stars decreases with increasing particle 
number. This trend shows that the high-mass stars, similarly to 
the finding for the density-scaled models, become less important 
as gravitational foci for the low-mass stars in clusters with larger 
stellar populations. 



7. Conclusion and Discussion 

The influence of different cluster environments on the encounter- 
induced disc-mass loss has been investigated by scaling the size, 
density and stellar number of the basic dynamical model of the 
ONC. The findings can be summarized as follows: 

1 . The disc-mass loss increases with cluster density but remains 
rather unaffected by the size of the stellar population. 

2. The density of the ONC itself marks a threshold: 

(a) in less dense and less massive clusters it is the massive 
stars that dominate the encounter-induced disc-mass loss 
via gravitational focusing of low-mass stars whereas 

(b) in denser and more massive clusters the interactions of 
low-mass stars with equal-mass perturbers play the ma- 
jor role for the removal of disc mass. 

3. In clusters four times sparser than the ONC the effect of en- 
counters is still apparent. 

These findings from numerical simulations are well con- 
firmed via observations. It is widely known that the disc fre- 
quency of even very young clusters (e.g. NGC 2024) is signifi- 
cantly below 100 % ( |Haisch et al.l200T]|Hillenbrand|2005| l. This 
implies a very short time-scale for the disc destruction of a part 
of the cluster population. In fact, stellar interactions are a po- 
tential candidate for a very rapid physical process and thus the 
encounter-induced disc-mass loss is potentially a vital mecha- 
nism for disc destruction at the earliest stages of cluster evo- 
lution. Observational evidence for the role of this mechanism 
has also been provided for the ~1 Myr old Orion Nebula Cluster 
( jOlczak et aTJ|20"08] l. Moreover, observations of NGC 2024, 
ONC, and NGC 3603 confirm also the correlation of decreas- 
ing cluster disc fraction (CDF) with increasing cluster density 
and even provide evidence for a critical density p cr j t ~ Ponc- 
A compilation of observational data shows that these similar 
aged (~1 Myr) clusters with peak densities of roughly 10 3 pc~ 3 , 



C. Olczak et al.: Stellar interactions in dense and sparse star clusters 



9 



family of models 






density-scaled 






size-scaled 


cluster model 


DO 


Dl 


D2 D3 


D4 


D5 


S0-S5 


low mass 


1/5 


2/5 


1 3 


7 


21 


1 


intermediate mass 


1/2 


1 


2 4 


9 


23 


2 


high mass 


5 


7 


10 15 


25 


44 


10 



Table 3. Approximate relative encounters rates f t 
low-mass group of the model D2/S2. 
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Fig. 9. Normalised encounter rate / e "°™ from the simulations in comparison with the analytical estimate given by Eq. \\5) . The 
filled squares represent all stars in the Trapezium Cluster region (R = 0.3 pc), the other open symbols stand for predefined mass 
groups: high-mass, m > 10 M Q (upward triangles), intermediate-mass, 10 M Q > m > 1 M Q (circles), and low-mass stars, m < 1 M Q 
(downward triangles). The lines depict the analytical estimate of the encounter rate for high-mass (dashed line) and low-mass stars 
(solid line). The ranges of the mass groups have been chosen here different from those in previous figures to account for the mass 
regimes of the encounter rate presented in Table 2] Left: Density-scaled cluster models. Right: Size-scaled cluster models (here the 
analytical estimates of the encounter rate for high-mass and low-mass stars are identical). 
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Fig. 10. Cluster density as a function of cluster size for clus- 
ters more massive than 10 3 M and embedded clusters with more 
than 200 observed members from Pfalzner (2009). The parame- 
ter space covered by the simulations in this work is indicated by 
the large pink cross centered on the ONC. 



10 4 pc~ 3 , a nd 10 s pc" 3 , respectively, show CDFs o f 85 %, 80 % 
and 4 % ( |Haisch et aLl|2000| |Lada et al.||2000| |Stolte et al.| 
[2004] i. The fact that NGC 2024 is sparser than the ONC but has 
a similar CDF is in good agreement with our simulations. 

Our results have several important implications for the gen- 
eral picture of star and cluster formation. Very recently, Pf alzner| 



(2009) has shown that there exist two cluster sequences evolv- 
ing in time along pre-defined tracks in the density-radius plane, 
the "leaky" and the "star-burst" clusters. The simulations per- 
formed in the present investigation cover the parameter space 
of the "leaky" clusters in their embedded stage (see Fig. 10 1. 



Comparison with our findings shows that at the earliest evo- 
lutionary stage leaky clusters have densities above the critical 
density. Hence in leaky clusters star-disc systems are initially 
efficiently destroyed via encounters that occur preferentially be- 
tween low-mass stars. The ONC corresponds to an intermediate 
stage in the embedded phase of leaky clusters with the transi- 
tion towards a preferred disc-mass loss of high-mass stars via 
gravitational focusing of low-mass stars. At the final stage of the 
embedded phase the encounter-induced disc-mass loss in leaky 
clusters ceases. Gravitational focusing by high-mass stars may 
still affect single objects yet at this age is most probably ex- 
ceeded by other disc-destructive processes like photoevaporation 
or planet formation. The effects in star-burst clusters would be 
similar yet much more pronounced. In case of the Arches clus- 
ter one could expect stellar encounters to destroy the discs of 
most of the low- and high-mass stars in several hundred thou- 
sand years. Combining our results with the finding of |Lada & 
|Lada| ( |2003| l that most stars are bom in clusters, it becomes ev- 
ident that a significant fraction of all stars must have been af- 
fected by stellar encounters at the early evolutionary stages of 
their hosting environment. 

The application of our results to the dynamics of embedded 
clusters - though obtained from simulations that do not contain 
gas - is justified for three reasons: 



10 



C. Olczak et al.: Stellar interactions in dense and sparse star clusters 



1 . Rather than simulating the evolution of leaky clusters (which 
would explicitly require the treatment of gas), we use our 
cluster models to map certain evolutionary stages of the se- 
quence of leaky clusters in terms of "dynamical snapshots". 
The dynamical effects of these numerical models are then 
used to estimate the effect of encounters in the observed clus- 
ters at their current dynamical state. 

2. The effect of gas in an embedded cluster is to lower the fre- 
quency of close encounters (due to the smoother cluster po- 
tential), yet unless the gas mass is dominating cluster dy- 
namics - as is not the case for the only partially embedded 
leaky clusters shown in Fig. 10 (e.g. NGC 2024) - the effect 
is minor. 

3. Gas expulsion causes the clusters to expand and thus their 
density to decrease much faster than in our simulations. 
Consequently, when mapping our results to the current dy- 
namical state, we underestimate the initial cluster density 
and thus the effect of encounters in the early evolutionary 
phases. Hence, our results are least accurate for older clus- 
ters, yet well applicable for the young evolutionary stages 
were encounters have the largest effect. 

The effect of the encounter-induced disc-mass loss in the 
early evolution of stellar systems has important implications for 
the formation and evolution of planets. From our findings we 
conclude that the probability to find planets around intermediate 
mass stars is highest. While the high initial densities of leaky 
clusters imply that planets around low-mass stars are expected 
to be less frequent, hardly any are expected in star-burst clusters. 
Independent of the cluster environment, planet formation around 
high-mass stars seems to be completely hindered via interactions 
with other cluster members. 
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Appendix A: Determination of boundaries of mass 
groups 

Boundaries of mass groups of low-, intermediate- and high-mass 
stars have been determined individually for different sizes of 
stellar populations on the basis of the IMF of Kroupa (2001 , see 
also Eq.[T]l. The derivation involves the requirement for the three 
mass ranges to be equidistant in logarithmic space, weighted by 
the slope of the IMF (of each mass range). The weighting ac- 
counts for the steepness of the slope in the high-mass regime 
which would otherwise cause a very sparsely populated group 
of high-mass stars. 
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In the case of a lower mass cutoff at mo = 0.08 M Q , and an 
upper mass limit m 3 , the IMF is characterised by just two dif- 
ferent slopes, ct\ = 1.3 in the range mo < m < 0.50 M , and 
ar2 = 2.3 in the range 0.50 M < m < m 3 . Because the break 
in the slope of the IMF at the critical mass = 0.5 M does 
not necessarily coincide with one of the boundaries of the mass 
ranges, the cases ni\ < m°J and m\ > m|f have to be differ- 
entiated. Though from the theoretical point of view the same 
differentiation would be required for the higher mass boundary 
TO2, this is not relevant for the stellar systems in the focus of the 
present work. The four mass ranges, m^, k — 0, .., 3, and the two 
slopes, at,k= 1,2, are then interrelated as follows: 



(m, > m°J) A (m 2 > m h J) : 

(log m\ - log m°J)a 2 1 + (log m h J - log mo)a\ 



1 



= (logw 2 - logmi)a 2 
= (logm 3 - logm 2 )Q'2 1 , 

{mi < m°J) A (m 2 > m]f) : 
(log mi - logmo)^ 1 
= (log m.2 - log m°J)ct2 1 + (log mjf - log m\)a\ x 
= (logm 3 - logm 2 )a2 1 . 

Solving these equations, and substituting 

a n = a\a 2 l , 
a 2l = u 2 a\ l , 
rn = 1 - a n , 
r 2 i = 1 - a 2 i, 

one obtains 



(m x > m^) A (m 2 > m|f) : 



log mi = 5 [logm 3 + 2r 2 i logm^ + 2a 2 i logm ] , 
log m 2 = 5 [2 log m 3 + T 2 i log m h J + a 2 \ log mo] , 



(A.l) 



(mi < m h J) A (m 2 > m h ") : 



log mi = 5 [ai2 log ra 3 + Fi2 log m°J + 2 log m J , 
logm 2 = 5 [21ogm 3 + T 2 i logmj! r + a 2 i logm ] . 

The choice of the appropriate solution is determined by the upper 
mass limit m 3 . For this purpose the "critical maximum mass" 
m™ x , 

< ax = log" 1 [(1 + 2an) log m h / - 2a l2 log m ] , 

is estimated from Eq. \A.l\ and m\ = m^ 1 . Consequently, the 
following relations hold: 

m-i < m™ x => mi > m\? , 

m 3 > w™ x => mi < m^ r . 

With the given values of the parameters mo, a\, and a 2 one finds 



The derived mass boundaries, m^, k = 0, .., 3, for each cluster 
of the families of models are presented in Table | A.l] 
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1000 


2000 


4000 


8000 


16000 


32000 


m [M o ] 


8.00 • 10~ 2 


8.00 • 10- 2 


8.00 • 10- 2 


8.00 • 10- 2 


8.00 • 10- 2 


8.00 • 10- 2 


mi [Mo] 


3.30-10-' 


3.54-10-' 


3.82- 10-' 


3.95 ■ 10-' 


4.15 ■ 10-' 


4.24 • 10-' 


m 2 [M ] 


2.94 


3.78 


4.95 


5.58 


6.61 


7.14 


m 3 [M ] 


1.47- 10 2 


1.47- 10 2 


1.48- 10 2 


1.48 • 10 2 


1.50- 10 2 


1.50 • 10 2 



Table A.l. Boundaries of the three mass groups of low-, intermediate-, and high-mass stars of cluster models with 1000, 2000, 
4000, 8000, 16000, and 32000 particles. 



